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OVERVIEW 

This  was  a  joint  project  with  MIT,  Henrik  Schmidt  et  al  (separate  ONR  proposal)  and  SACLANTCEN 
(Stevenson  and  Jensen).  The  long  term  goals  are  to  develop  environmentally  adaptive  bi-  and  multi¬ 
static  sonar  concepts  for  autonomous  off-board  sensor  networks  for  the  detection  and  classification  of 
proud,  buried  and  waterborne  targets  in  shallow  water.  SACLANTCEN  provided  three  weeks  of 
simultaneous  sea  time  for  both  the  R/V  Alliance  and  R/V  Leonardo  during  July  2004.  MPL’S  main 
contribution  was  to  develop  and  test  signal  processing  algorithms. 


OBJECTIVES 

The  objective  of  the  research  proposed  was  to  combine  the  results  of  two  ONR  sponsored 
SACLANTCEN  Joint  Research  Projects  (JRP),  the  GOATS  (Generic  Oceanographic  Array  Technology 
Sonar)  multistatic  MCM  initiative  led  by  MIT,  and  the  FAF  (Focused  Acoustic  Fields)  time  reversal 
initiative,  led  by  MPL,  into  a  new,  environmentally  adaptive,  multistatic  sonar  concept  for  concurrent 
detection  and  classification  of  proud,  buried  and  waterborne  targets  in  shallow  water,  using  off-board 
autonomous  sensing  platforms  in  conjunction  with  organic,  on-board  sonar  resources.  The  research  was 
to  develop  a  comprehensive  System  Modeling  and  Simulation  framework  which  will  be  used  for 
performance  prediction  and  experiment  planning.  Concept  demonstration  was  to  be  performed  through  a 
FY04  field  experiment  conducted  jointly  by  the  three  organizations. 
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SONAR  CONCEPT 


Figure  1.  Multistatic,  low-frequency  sonar  concept  for  concurrent  detection  and  classification  of  proud  and  buried 
targets  in  shallow  water,  using  off-board,  autonomous  multi-static  sonar  platforms  in  conjunction  with  organic 
on-board  sonar  resources  and  time-reversal  focusing. 

The  sonar  concept,  shown  schematically  in  Figure  1,  is  derived  from  GOATS,  the  multi-static, 
autonomous  MCM  concept  for  very  shallow  water  currently  being  explored  by  MIT  under  ONR  funding 
(Codes  321  and  322)  in  collaboration  with  SACLANTCEN.  It  uses  a  fleet  of  autonomous  off-board 
platforms  to  create  a  wide,  combined  physical  and  synthetic  aperture  for  mapping  the  spatial  and 
temporal  characteristics  of  the  scattered  field  produced  by  the  targets  in  the  mid- frequency  regime 
which  is  rich  in  target-specific  structural  responses.  The  multistatic  field  characteristics,  rather  than 
classical  sonar  imaging,  is  then  used  as  the  the  basis  for  concurrent  detection  and  classification.  For 
insonification  the  concept  would  be  to  use  organic,  mid-frequency  sonar  resources  on  a  littoral  surface 
platform,  e.g.  the  SQS53  or  a  wider-band  successor  on  the  Littoral  Combat  Ship  (LCS),  in  an 
environmentally  adaptive,  time-reversal  mode.  Thus,  the  time-windowed  reverberation  is  used  to 
determine  the  instantaneous  impulse  responses  from  a  several  kilometer  wide,  broadside  swath  of 
seabed,  capturing  all  the  shallow  water  multipath  structure  under  the  current  environmental  conditions. 
These  inpulse  responses  are  then  convolved  with  the  desired  waveform,  time-reversed  and  transmitted 
by  the  sonar  to  produce  a  focused  field  on  the  seabed,  with  the  focus  scanned  over  the  entire  swath,  out 
to  ranges  of  several  kilometers.  Algorithm  development  for  such  use  of  time  reversal  reverberation  for 
enhanced  sonar  performance  is  currently  being  explored  under  ONR  funding  (Code  321  US,  Doug 
Abraham). 

A  network  of  autonomous  underwater  vehicles  are  used  to  provide  an  off-board,  multistatic,  synthetic 
receiving  array  and  bi-static  SAS  concepts,  developed  under  GOATS,  to  detect,  track  and  classify  proud 
and  buried  targets  in  the  insonified  sonar  footprint.  Once  detection  is  established,  the  target  coordinates 
are  transmitted  back  to  the  surface  platform,  which  will  subsequently  focus  the  sound  on  the  target 
location,  while  the  bistatic  platforms  adaptively  change  their  path  to  optimally  classify  the  target  from 
the  spatial  and  temporal  characteristics  of  the  bistatic  scattering,  using  adaptive  path-planning 
techniques  developed  and  demonstrated  under  GOATS.  The  AUV’s  also  carry  broad-band  sources 
which  may  be  applied  for  close-in  target  insonification,  using  adaptive  path-planning  and  possibly  time- 
reversed  target  replicas  for  rapid  classification. 
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APPROACH 


MPL’s  role  in  this  program  concerned  signal  processing  for  target  detection.  Hence  we  were  to  develop 
and  demonstrate  a  self-adaptive  focusing  procedure  that  required  minimial  enviommental  information  in 
order  to  scan  the  bottom  for  targets.  This  was  first  explored  in  the  laboratory  and  then  subsequently  at 
sea. 

RESULTS 

Self-adaptive  Focusing. 

We  developed  an  iterative  time  reversal  procedure  to  focus  on  a  target/  The  first  experiments  were  done 
in  a  laboratory.  After  demonstrating  the  procedure  in  the  laboratory,  data  was  taken  at  sea  and  analysed 
indicating  that  the  procedure  does  indeed  work.  The  laboratory  setup  is  shown  in  Figure  2. 


b) 


•Ultrasonics  linear  array, 

•Frequency  Bandwidth  [1Mhz-2Mhz] 

•Ic  =1mm  (at  Fc=1.5Mhz, ) 

•M=12  elements,  equally  spaced  by  2.5l 

•Elements  size  (12mm*0.5mm)->high 
angular  directivity-  >2D  acoustic  progation 
mainly 


Figure  2.  Laboratory  setup  for  ultrasonics  experiment  to  demonstrate  iterative  time  reversal. 

The  iterative  time  reversal  method  shown  schematically  in  Figure  3  was  the  basis  of  what  was  used  to 
construct  the  reflection  maps.  However,  during  practical  mapping  of  the  ocean  bottom,  the  TRA  is 
moving  between  each  iteration  of  the  active  time-reversal  process  since  it  is  either  being  towed  by  a  ship 
or  directly  mounted  on  the  ship’s  hull.  This  TRA  motion  is  detrimental  since  the  TRA  would  not 
retrofocus  on  the  exact  same  range  section  on  the  bottom.  We  use  a  passive  implementation  of  the 
iterative  time-reversal  process  to  overcome  this  issue.  This  requires  the  knowledge  of  the  interelement 
impulse  response  matrix  Ki,mof  the  M=12  elements  of  TRA:  each  element  ki>m(t)  of  the  matrix  Ki>mis  the 
output  of  the  element  number  1  (1<1<M),  when  the  input  element  number  m  (l<m<M)  is  a  delta  impulse 
[1,3,4],  In  other  words,  the  kiM(t)  signals  correspond  to  the  backscattered  Green’s  functions  between  the 
1th  and  mth  element  of  the  TRA,  including  the  bottom  reverberation,  convolved  with  the  probing  pulse 
signal.  Then,  the  time  reversal  iterations  are  implemented  passively  by  simply  performing  correlations  of 
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these  k  i>m(t)  signals  with  the  time-gated  window  of  the  reverberated  signals  selected  at  the  previous 
iterations. 


Emission  1  Reception  1 

STEP1 


broadside  illumination 


transducer  array 


_  srattprpr 


backscattered  signals 


STEP  2 


time  window  selection 


^  * jj 

time  reversal 


time  reversed  window 


‘  ’  -  ft* 

£ 

■JIT'' 

STEP  3 


Emission  2 

time  reversed  window 


backscattered  signals 


Reception  2 


Figure  3.  Schematic  of  iterative  time  reversal 


The  interelement  impulse  response  matrix  Klmcan  be  initially  measured  by  firing  successively  the  same 
short  impulse  signal  from  each  individual  array  element  and  then  recording  the  backscattered  signals  on 
the  whole  TRA  after  a  round-trip  propagation  between  the  TRA  and  the  ocean  bottom.  To  avoid  waiting 
for  M  succesive  reverberation  returns  to  be  completed  (which  might  be  too  long  if  the  TRA  is  moving), 
it  is  possible  to  record  in  one  shot  the  matrix  K,  m  by  probing  the  waveguide  with  a  set  of  M  orthogonal 
array  beams  rather  than  the  individual  TRA  elements. 

For  each  of  the  kiM(t)  signals,  the  recording  time  t  corresponds  to  the  round  trip  propagation  from  the 
TRA  to  the  water-bottom  and  can  thus  be  converted  to  an  equivalent  distance  d=tco  (in  wavelengths,  X), 
where  Co  is  the  speed  of  sound  in  the  water  tank.  The  reflectivity  maps  along  the  waveguide  bottom 
were  constructed  using  the  following  steps  (see  Fig.4): 

1 .  Measure  the  interelement  impulse  response  Ki>m  of  the  TRA,  using  a  chirp  signal.  This  is  the 
only  active  part  of  the  process. 

2.  From  Ki,m  ,  build  passively  the  M  bakcscattered  signals  R(0\d)  (l<m<M)  issued  form  a 
broadside  transmission  (see  Fig.  3) 
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(1) 


M 

K V)=IX(<o 

i=\ 

3.  Select  several  time-gated  windows  (of  spatial  length  Aco),  each  associated  with  a  different 
average  range  on  the  bottom  Lav  determined  by  the  center  of  selected  time-gated  window  (see 
Fig.  3,  for  Lav=105A.  where  the  window  length  zl=5ps  is  delimited  by  black  lines) 

4.  Perform  passively  the  iterative  time-reversal  process  for  each  window,  by  simple  correlations 
(indicated  by  the  symbol  0 )  with  the  elements  of  Ki>m,  for  a  chosen  number  of  iterations  N 
(see  Fig.  3.b). 

M  A  r  A  r 

R'‘>(d)  =  £*„(</)« <  d  < 4,  +— i)  (1  <  i  <  N) (2) 

/= i  4  1 

5.  Plot  the  average  recorded  energy  by  the  TRA,  Eav(d,Lav),  vs.  the  illuminated  average  bottom 
range  Lav  (see  Fig.  1,  Fig.  5  and  Fig.6).  Eav(d,Lav)  is  computed  incoherently,  by  simply  time¬ 
aligning  the  maxima  of  the  envelope  of  each  retrofocused  energy  signals  R(N\d )  before 
doing  an  incoherent  summation. 


Step  1:  Recording  of  the  inter-element  response  matrixK|  m 
(use  set  orthogonal  signal^Mast  acquisition) 


Elt  #1  broadcast  l(t) 
Elt  #m  records  k  lm(t) 


k  im(d)  =  backscattered  Green ’s  function  between  the  Ith  and  mth  element  of  the  TRA 
including  the  bottom  reverberation  ,  convolved  with  the  probing  pulse  signal. 

Step  2:  Passive  time  reversal  iterations:  use  correlations  of  thesek 
l,m(t)  signals  with  the  time-gated  window  of  the  reverberated  signals 
computed  at  the  previous  iterations  R(° 

M  El  El 

*1V)  =  --<d<L„+-)  (1  <i<N) 

1=1  4  2 


Figure  4.  Simplified  schematic  of  the  “passive”  implementation  of  iterative  time  reversal. 

An  example  of  laboratory  results  for  a  reflectivty  map  with  a  target  is  shown  in  Figure  5  where  four 
iterations  were  carried  out. 
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Envelope  of  the  average  recorded  energy  @  the  TRA 
EJd,Lav=W5/) 

Figure  5  Ultrasonic  Laboratory  results  indicating  target  is  acquired  by  a  self-adaptive  process. 

At-Sea  experimental  results 

Among  the  arrays  deployed  in  the  FAF04  experiment  was  the  “billboard”  arrays  shown  in  Figure.  6  . 


Figure  6:  Planar  array  used  in  FAF04 
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An  example  of  results  obtained  with  the  Billoboard  array  are  shown  in  Figure  7  in  which  there  is  also  a 
comparision  with  a  noniterative  result  and  a  broadside  non  adaptive  result.  Clearly,  the  iterative  TR 
result  inhances  the  target  in  the  reflectivity  map. 


TARGET  in  the  Billboard  axis 
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Figure  7.  Experimental  results  for  an  at-sea  reflectivty  map  showing  target  detection. 
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